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We study quantum teleportation via Bell-diagonal mixed states of two qubits in the context of 
the intrinsic properties of the quantum discord. We show that when the quantum-correlated state 
of the two qubits is used for quantum teleportation the character of the teleportation efficiency 
changes substantially depending on the Bell-diagonal-state parameters, which can be seen when the 
worst-case-scenario or best-case-scenario fidelity is studied. Depending on the parameter range, one 
of two types of single qubit states is hardest/easiest to teleport. The transition between these two 
parameter ranges coincides exactly with the transition between the range of classical correlation 
decay and quantum correlation decay characteristic for the evolution of the quantum discord. The 
correspondence provides a physical interpretation for the prominent feature of the decay of the 
quantum discord. 


I. INTRODUCTION 


The quantum discord , a measure of bi- and multi¬ 
partite quantum correlations, has attracted much at¬ 
tention recently. This is due to the fact that the dis¬ 
cord indicates the existence of quantum correlations in 
many partially mixed states which have no entanglement 
present. Specifically, the quantum discord, V , does not 
display any sudden death type phenomenon, since the 
zero-discord states form a set of measure zero fl]. This 
subset of states is the set of all truly classical states 
chosen asymptotically by decoherence (SI, 0. Hence, a 
smooth, continuous decoherence process cannot lead to 
a sudden and continued disappearance of quantum cor¬ 
relations mid-evolution (before a fully mixed, completely 
dephased state is reached). This suggests that the sud¬ 
den death of entanglement signifies not, as was previously 
believed, the disappearance of all quantum correlations, 
but the crossing of a threshold of a given, small amount 
of correlations and the disappearance of quantum corre¬ 
lations of a certain type. Although below this threshold 
many quantum informational tasks are no longer pos¬ 
sible, methods of performing quantum computation on 
zero-entanglement states for which the quantum discord 
is non-zero have already been devised |7H11|. Further¬ 
more, it has been very recently shown experimentally 
how entanglement can be shared between distant parties 
via non-entangled states with non-zero discord fl2| - (l6j . 


Unfortunately, computing the quantum discord for an 
arbitrary density matrix is an extremely involved task 
even in the simplest two-qubit case 17]. This led to the 
emergence of the geometric quantum discord jl8| , which 
is defined as the minimal Hilbert-Schmidt distance of a 
given state from the set of zero-discord states. Although 
an explicit formula for the geometric discord given a two- 
qubit density matrix does not yet exist (one that does 
not require minimization over the set of all zero-discord 
states), such formulas exist for the lower [l8| and up¬ 
per [191 ] bounds on the geometric discord. The geometric 


discord is a good measure to distinguish between zero- 
discord and non-zero-discord states, but because of the 
properties of the Hilbert-Schmidt distance, it is not a 
good measure for the amount of quantum correlations 
present in a given state. In fact, because the Hilbert- 
Schmidt distance is sensitive to the purity of the state, 
the geometric discord may be increased by non-unitary 
evolution of a single qubit (the unmeasured one) [201 - 122 , 
which should not increase inter-qubit quantum correla¬ 
tions. One solution to this problem is rescaling the ge¬ 
ometric discord following Ref. (23}; a lot of work is also 
being done to find explicit formulas for discord measures 
utilizing more appropriate distance measures [24l - l27i ]. 


Regardless of these difficulties, there are properties of 
typical discord evolutions, T>(t), under decoherence pro¬ 
cesses that have already become apparent. Perhaps the 
most prominent and most baffling is the transition be¬ 
tween the regimes in which either the classical correla¬ 
tions, C (defined as the difference between the mutual 
information and the quantum discord T>), or the quan¬ 
tum correlations (i.e. V itself), exhibit a faster decay pgH - 
[30| . The most striking example of such a behavior was 
given in [29 ], where it was shown that for a certain ini¬ 
tial Bell-diagonal state, the influence of phase-damping 
channel leads to the evolution, during which T> remains 
constant for a finite time t c after the state initialization 
(while C(t) are decaying), and for longer times it decays 
(while the C{t) remain constant). In a more general case, 
the transition between the two regimes is manifested by 
discontinuity in the first derivative of T>(t ), i.e. by a tran¬ 
sition between two regimes governed by different decay 
functions. Such a behavior was also seen in the dynamics 
of the geometrical quantum discord and rescaled geomet¬ 
ric quantum discord [23], and the transition was shown 
to occur at exactly the same point as the transition in the 
decay of V(t) (31, [Hj. The study of the geometric dis¬ 
cord revealed similar transitions occurring for non-Bell- 
diagonal initial states and evolutions. 


The above discussion suggests that the transition be- 




2 


tween the “classical” and “quantum” decay regimes is 
a rather general feature of T> dynamics. There is, how¬ 
ever, no intuitive understanding of the possible opera¬ 
tional meaning of this transition, which is not surprising 
in the light of the fact that the operational significance 
of the quantum discord is a subject of ongoing research 
[33l [34l ]. In the following, we give a simple interpreta¬ 
tion of an observable (and practically relevant) quantity 
that changes discontinuously at the transition point of 
discord dynamics. We investigate the quantum telepor¬ 
tation protocol with the entangled Bell state used for 
teleportation being subject to decoherence. It is well 
known that the time at which the average teleportation 
fidelity becomes smaller than 2/3 is the time at which the 

t artite state used as a resource becomes separable (2^- 
. Here we focus on teleportation fidelity minimized 
over the teleported states (i.e. the worst-case-scenario fi¬ 
delity) or maximized over the teleported states (i.e. the 
best-case-scenario fidelity; the two situations display an 
equivalent transition), and we show that at the time of 
the transition between the two regimes of quantum dis¬ 
cord decay, the nature of the state for which teleportation 
has the lowest/highest fidelity changes. 


where {n^} is a complete set of orthonormal projective 
operators corresponding to a von Neumann measurement 
of subsystem B. The index k denotes the outcome of a 
given measurement and the formula involves maximiza¬ 
tion over the set of projective measurements. There¬ 
fore the formula of Eq. (© yields the information gained 
about the system A after the measurement {Hfc} on sys¬ 
tem B. The quantum discord of a given state pab is then 
given by 


T^{pab) = I{pab) ~ C(pab)- (3) 

Following Ref. [HI, the value of the quantum discord 
can be found for any Bell-diagonal two-qubit state. If the 
Bell-diagonal two-qubit density matrix is written in the 
form 


PAB = 


( P<s> 
0 
0 

V 0-$ 
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where all four parameters are real, 
tion is given by 
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II. CLASSICAL AND QUANTUM 

DECOHERENCE REGIMES OF QUANTUM 
DISCORD DYNAMICS 

The quantum discord is defined as the difference be¬ 
tween two classically equivalent formulas for mutual in¬ 
formation [l] . The formula which is referred to as mutual 
information in Ref. [29] is given by [lj 

Z(pab) = S(Tt a Pab) + S(Ti b pab ) - S(pab), (1) 

where the von Neumann entropy is given by S(p) = 
— Trplog 2 p. This quantity was generalized from 
the classical language of probability distributions in a 
straightforward manner to the language of density ma¬ 
trices, while the Shannon entropy was replaced by von 
Neumann entropy. The other formula for classical mu¬ 
tual information, which is in the quantum context re¬ 
ferred to as classical correlations (in Ref. (29j), cannot 
be generalized in a direct manner, because the classical 
formula involves conditional entropy, 


1(PAB) = 2+ ^ l(Pi + 0i)l0g 2 (pi + CTi) 

+{Pi - d) log 2 {pi - at)}. (5) 

Note that any continuous and differentiable evolution 
of pab that retains the Bell-diagonal form of Eq. © 
must lead to a continuous and differentiable evolution 
of I(pab)- An explicit formula for the classical correla¬ 
tions of Eq. © can also be found for any Bell-diagonal 
state [HI . These correlations are given by 

C(pab) = ^ K 1 + X) log 2 (l + X) + (1 - X) log 2 (l - X)], 

( 6 ) 

where x = max{|A|, \<Jm\ + |cr<j>|}, in which we used 
A = pm — pm- The maximization allows for indifferen¬ 
tiability points in the evolution of classical correlations 
which occurs when the plane 

|A| = |cNr| + |f7$| (7) 


C(A : B) = H(A) — H(A\B), 

where H denotes the Shannon entropy and A and B are 
random variables. Conditional entropy H(A\B) requires 
the specification of the state of A given the state of B, 
which in quantum mechanics is ambiguous until the mea¬ 
surement performed on B is specified. Hence, the con¬ 
ditional von Neumann entropy can be found given the 
complete measurement on subsystem B and the result¬ 
ing formula for classical correlations is [h, 0 

C(pab) = max [S(Tr pab) ~ S/pAsKlIfc})], (2) 

{n fc > 


is transgressed, resulting in indifferentiability points of 
the quantum discord which is given by the difference 
of the smooth mutual information and the indifferen- 
tiable classical correlation function, Eq. ©■ Note that 
under a pure dephasing decoherence process, the clas¬ 
sical correlation function remains constant in the |A| > 

|< 7 ^ | + |< 7 $ | regime, which we will call the quantum de¬ 
coherence regime following Ref. [HI, while it decays in 
the classical decoherence regime, |A| < |cr$| + |cr$|. Be¬ 
low we consider decoherence which leads to C(t) and T>(t) 
both decaying in the two regimes defined by the above 
inequalities. 
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III. RELATION BETWEEN TELEPORTATION 
FIDELITY AND DISCORD 

To gain some understanding of the physical significance 
of this transition, let us turn to the quantum teleporta¬ 
tion of an unknown qubit state by means of an entangled 
two-qubit state [38lj in the situation when the entangled 
state has previously undergone partial decoherence. We 
will focus on the scenario where the initial entangled state 
is the |<f> + } = l/\/2(|00) + 111)) Bell state, corresponding 
to = cr$ = 1/2 and p-q, = crq, = 0 in Eq. (|5j); the as¬ 
sumption is only made for simplicity; the same results are 
acquired regardless of the chosen Bell state. For such an 
initial state, decoherence often retains the Bell-diagonal 
form, but rarely provides the means to induce coherences 
between the |01) and 110) components (cr^^O). It can, 
however, induce non-zero occupations corresponding to 
|01) and 110) [SflL l4(ij|. and possibly change the diagonal 
matrix elements in such a way that the state ceases to 
be Bell-diagonal. Since we want to maintain the Bell- 
diagonal form at all times (so that we can analytically 
calculate the quantum discord), below we will focus on a 
physically well motivated example of decohering channel 
which preserves this form with <jq,=0. 

First let us explain the connection between the two 
regimes of discord dynamics and the teleportation fi¬ 
delity. The unknown qubit to be teleported is \ip) = 
a|0) + /3|1), with \a\ 2 + \(3\ 2 = 1. After teleporting | ip) 
using the state of Eq. with cr^i = 0 instead of the 
|<h + ), the teleported state is equal to 

Pt = |0)(0|(2p$|a| 2 + 2p^|/3| 2 ) (8) 

+ |1}(1|(2 P *H 2 + 2p*|/3| 2 ) + |0><l|2<7*a/3* + H.c.. , 

and the fidelity of teleportation is given by 

F 2 = {$\ptW) = 2- 4(A - cr$)|a| 2 (l - |a| 2 ) . (9) 

The standard procedure now is to find the average fidelity 
of the teleported state (averaged over all possible states 
to be teleported) and compare it to the classical limit 
of teleportation capability which is equal to 2/3 (35l-l37i|. 
The average fidelity is given by 

F 2 V = ^P<s> ~ ^(A — cr$) , (10) 

which is a smooth function of time. In order to see the 
correspondence between the transition in the evolution 
of quantum correlations and teleportation capability, one 
should rather look at the extremal scenarios, namely, the 
fidelity of the teleported state minimized or maximized 
over all the | i/j) states. Analyzing the extrema of Eq. <[9|) 
as a function of |a| 2 G [0,1] we find that for |A| < |cr$| 
(in the classical decoherence regime), the minimal fidelity 
^min = 1/2 +A occurs for | ip) = |0), |1) (the poles on the 
Bloch sphere), and the maximal fidelity F^ ax = l/2 + er$ 
is obtained for \ip) = l/\/2(|0) + exp(*(/>)|l) (the states 
on the equator of the Bloch sphere). On the other hand, 
for |A| > |(7$| (in the quantum decoherence regime), the 


states which are easiest/harder to teleport (with the cor¬ 
responding fidelities) trade places, so that we always have 
+A n = l/2 + min(A,er$) and ax = 1/2 + max(A, cr$). 
At the transition point, A = cr$, the fidelity from Eq. © 
is independent of the state \il>), so that it is equal to its 
average value, 2p$ = A + 1/2. If A >1/6 at this point in 
time, the transition between the two regimes of quantum 
discord dynamics occurs when the state is still entangled, 
and F 2 V is larger than its maximal classical value. 

The character of state \(f>) which is easiest/hardest to 
teleport in a given regime can be also connected with the 
character of the classical states closest to pab EH. In 
the “classical” regime there are two such states, p x and 
p v given by 

Pa = (l + t) (MM + MM) + 

(l ~ t) (MM + MM) (u) 

where a = x,y , and the states |a) and |a) are eigenstates 
of the Pauli matrices <j a with ±1 eigenvalues. Such a 
classically correlated state p a allows for “teleportation” 
of \(j>) = |a), |a) with fidelity of 1/2 + cr$. The fact that 
these are the closest classical states in cr$ > A regime 
gives an intuitive explanation for the maximal fidelity 
being achieved for \cf>) state having its Bloch vector in the 
xy plane. On the other hand, in the “quantum” regime 
the closest classical state is given by formula analogous 
to the one from Eq. m , but with a = z and cr$ being 
replaced by A. The correlations in this state allow for 
classical enhancement of teleportation fidelity of |0) and 
|1) states. 

Fig. |T] shows the regimes of decoherence for Bell- 
diagonal density matrices with aq, = 0 with respect to 
the amplitude of the coherence present in the system and 
the difference of the two distinct two-qubit occupations. 
The border between the region of quantum decoherence 
(cyan) and classical decoherence (yellow) in terms of the 
quantum discord is denoted by the blue line. The Bell 
state used above is located in the upper right corner of 
the figure (A = 1/2 and |cr$| = 1/2) and its decoherence 
will move it to the lower left. In fact, a pure dephas¬ 
ing process of an initial Bell state results in its moving 
straight down along the A = 1/2 line, which is located 
solely in the quantum decoherence regime. Vanishing 
entanglement corresponds to the point when F 2 reaches 
the classical communication limit of 2/3 )35li371 ]. This 
boundary is denoted by a red, dashed line in Fig. [Q As 
seen, sudden death of entanglement EH"EI1 (i.e. crossing 
of the separability boundary while (7$ is nonzero) can 
be achieved both via the quantum decoherence regime 
and the classical decoherence regime. Moreover, entan¬ 
glement sudden death can occur for higher values of co¬ 
herence cr$ in the “classical” parameter range, because 
it is susceptible to the disturbance of qubit occupations 
EH- Incidentally, the zero-discord line is the |cr$| = 0 
axis and coincides with the minimized teleportation fi¬ 
delity reaching the minimal value of 1/2. 
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FIG. 1. Regimes of decoherence depending on the Bell- 
diagonal density matrix parameters for < 7 * = 0. The blue 
line denotes the indifferentiability plane given by Eq. (f7|l . di¬ 
viding the quantum decoherence regime (cyan triangle) from 
the classical decoherence regime (yellow triangle). The un¬ 
colored region corresponds to parameter values that yield 
an unphysical density matrix. The red dashed line is the 
border, zero-entanglement line (entanglement vanishes below 
the line). The dotted lines are izodiscords, corresponding to 
D = 0.1, 0.3, 0.5, 0.7 going from bottom to top. 

It is straightforward to generalize the teleportation 
procedure via a decohered two-qubit state to any Bell- 
diagonal state, Eq. ©. and to show that the discord 
transition plane of Eq. 0 still corresponds to a tran¬ 
sition in the worst-case-scenario teleportation from the 
situation when the states on the pole of the Bloch sphere 
are hardest to teleport, and when states on the equator 
of the Bloch sphere are hardest to teleport. Hence the re¬ 
lation between quantum teleportation and the quantum 
discord presented in this Section holds for any state of 
the form given by Eq. ©• 

IV. POSSIBLE PHYSICAL REALIZATION 

Let us now present a model of open-system dynamics 
in which the above-discussed transition between two tele¬ 
portation regimes occurs at finite time t c . It is clear that 
the classical decoherence requires a minimal amount of 
dephasing (decay of <r$) consistent with a given amount 
of relaxation (decay of A), suggesting the use of am¬ 
plitude damping channel. In order for the state to re¬ 
main Bell-diagonal, we use the generalized amplitude 
damping channel [39| corresponding to an environment 
being in equilibrium at temperature much higher than 
the energy splitting of the qubits. Assuming Marko¬ 
vian decoherence, we have then A(t) = A(0)e~ rt and 
o$(f) =cr$(0)e -rt / 2 , where r = TA + Tb, with T a,b be¬ 
ing the longitudinal relaxation rates of the two qubits 
forming the entangled state used for teleportation. 

However, under the action of this channel alone, the 


evolution starting form one of the Bell states remains in 
the classical regime at all times. In order for the tran¬ 
sition to the quantum decoherence regime to occur, the 
additional dephasing process, having negligible contribu¬ 
tion at short times, has to become stronger at longer 
times. This will happen when the qubits forming the 
entangled state used for teleportation are coupled via 
their d z operators to a slowly fluctuating bath. Low 
frequency (w <C where T is the bath tempera¬ 

ture) environmental fluctuations can be treated as clas¬ 
sical noise £(t) [HI, which under a realistic assumption 
of Gaussian statistics is described by a spectral density 
S(w) = /e lajt (£(£)£(0))df (where (...) denotes averaging 
over noise realizations. When most of spectral weight is 
concentrated at low frequencies (u> -C I/T 2 where T 2 is 
the characteristic dephasing timescale), the total noise 
power is a 2 oc Jq^ 2 S(u>)dcu. This is a typical situa¬ 
tion in many solid-state based qubits (with 1/f charge 
and flux noise 147 and hyperfine interaction with nu¬ 
clear bath [40j, l48l l49j being prominent examples). Such 
a noise leads to an additional decay of er$(f), given by 
a $(£) —>■ cr$(f) exp[— (q£) 2 ], with j = cr/y/2. Note that 
since the dephasing and relaxation are caused by very 
different processes (the former by low-energy fluctuations 
in the environment, the latter by high-energy processes 
involving exchanges of energy quanta corresponding to 
qubits’ energy splittings, which are typically ^ ksT), 
treating them as independent and additive processes is 
a reasonable approximation. Decay of quantum correla¬ 
tions due to such low-frequency or quasi-static baths was 
a subject of a few recent works [321. 150l - l54| . Within this 
model of decoherence we obtain a transition between two 
regimes of quantum discord decay at £ C = T /2q 2 , at which 
point we have F^ in = i[l + exp(— |(T/q) 2 ], so that this 
worst-case scenario fidelity still exhibits entanglement- 
related enhancement when T/q < y/2 In 3. The modulus 
of the discontinuity of the derivative of 1~>(t) at t = t c 
is given by ±TA c ln(^§^), where A c = |exp(-T/q). 
The quantum-classical transition is thus most visible for 
low values of T/q, which means that the rms a of the 
low-frequency phase noise should be much larger than 
the energy relaxation rate of the qubits. When the en¬ 
ergy relaxation is in fact caused by transverse coupling to 
high-frequency noise, the power of this noise at frequen¬ 
cies corresponding to qubits’ energy splittings (which is 
proportional to T) should be much smaller than a. When 
the low- and high-frequency noises have the same physi¬ 
cal origin, this requirement is consistent with our assump¬ 
tion of the low-frequency character of the environmental 
fluctuations. 

Note that while the above result has been obtained 
in the case of pure dephasing coupling to low-frequency 
noise (which leads to Gaussian decay of cr$(t)), a quali¬ 
tatively analogous transition between classical and quan¬ 
tum decoherence regimes can occur for more general cou¬ 
pling to low-frequency noise. For example, for transverse 
coupling to noise in the presence of large energy split¬ 
ting (i.e. the qubit Hamiltonian H = fla z /2 + £(t)a x /2 
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with fl 7 > a/(£ 2 )~), the evolution of the qubits is approx¬ 
imately [55h59I | of the pure dephasing form (only with 
the term quadratic in the noise, £ 2 (f)/ 2 f 2 , contributing 
to the qubits’ splittings), and the initial coherence decay 
is a$(t) ~ 1 — {t/"/') 2 , with 7 ' « 7 2 /Q. At very short 
times (t/ 7 ') 2 -CTf, and the discord dynamics will be in 
the classical regime, while at longer times the quantum 
decoherence regime can be entered. Interestingly, since 
the asymptotic decay of a$(t) in this case is of power- 
law kind 0 , Ull HI] (so that for long t again we have 
cr$(f) > A(t)), the transition into the quantum decoher¬ 
ence regime has to be followed by a re-entry into the 
classical decoherence regime at a later time. 

V. CONCLUSION 

Summarizing, we have studied the properties of quan¬ 
tum correlations present in mixed two-qubit states of 
Bell-diagonal form. We have shown that there is a direct 
correspondence between the transition point between the 
regime of classical correlation decay and the regime of 
quantum correlation decay, which is a characteristic of 
the evolution of the quantum discord, and the transi¬ 
tion point which appears for the worst-case-scenario Fi¬ 


delity of quantum teleportation. For teleportation, the 
two regions being transgressed correspond to two differ¬ 
ent classes of states which are hardest to teleport. Those 
are either the equal superposition states in the quan¬ 
tum decoherence regime, or the single-component | 0 )/| 1 ) 
states in the classical decoherence regime. This shows 
that there is a qualitative physical difference between 
the quantumly correlated states studied, depending on 
the correlation decay regime in which the state is lo¬ 
cated. Furthermore, the fact that the transition point 
is seen also by studying the geometric discord and that 
similar transition points are seen for initial states which 
for technical reasons cannot be handled using the regu¬ 
lar quantum discord, suggests that there is an underlying 
physical meaning to all such transition points, and that 
classical-quantum decoherence transition is a widely oc¬ 
curring property of the quantum discord. 
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